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| Abstract. We propose a method to investigate the scenario that cosmic relic 

neutrinos are highly clustered around stars and galaxies, or dark-matter clusters, 
rather than uniformly distributed in the universe. Such a scenario can be detected 
or constrained by the interaction of high energy cosmic ray protons and nuclei with 
the cosmic neutrinos. There should be observable signature in the energy spectra 
, of cosmic ray protons and nuclei for a neutrino clustering factor beyond 10 13 . We 

' provide a relation on the signature onset positions between proton and nuclei spectra, 

C*~) • and discuss possible support from existing experiments. It is also suggested that the 

relative abundance of cosmic ray nuclei may detect or constrain the cosmic neutrinos 
' with smaller clustering. 
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From the standard models of cosmology and particle physics, the history of our 
universe can be traced back to about 10~ n seconds after the Big Bang t L using well- 
established physical laws. At that time the universe was a mixture of photons, leptons, 
and quarks at thermal equilibrium with temperature of T « 300 GeV. Neutrinos 
(antineutrinos) decoupled from the thermal universe at the universe age of 1 second 
with T K 1 MeV, and photons decoupled from the baryonic matter at the age of 
3 x 10 5 years with T « 0.3 eV. Afterward, neutrinos, photons, and baryonic matter 
are supposed to evolve independently as the universe expands. The mean number 
densities of these three components in the present universe are well estimated as 
n v = (3/22)n 7 « 56 cm~ 3 per species of neutrinos (antineutrinos), n 7 = 413 cm -3 
for photons, and ng~5x 10~ 10 n 7 w 2 x 10~ 7 cm~ 3 for baryonic matter. From well 
established observations of cosmic microwave background (CMB) [2 , we know that 
photons are distributed in uniform density with small fluctuation of the order 10~ 5 |3j. 
However, baryonic matter are distributed in clusters as galaxies and stars, as a result of 
gravitational attraction. Taking the sun as an example, the baryon number density is 
n Q = 8.4 x 10 23 cm -3 , with a clustering factor of Uq/ub ~ 4 x 10 30 compared with the 
mean baryon number density of the universe. Thus, photons and baryonic matter have 
distinct patterns of distributions in the universe. 

It is completely unknown whether neutrinos are distributed in uniform density 
in our universe or gathered as clusters around the stars and galaxies. One would 
expect a uniform distribution in analogy to photons, if the neutrino mass is exact 
zero. However, it has become evident that neutrinos do have nonzero masses from 
the oscillation measurements jlj. So one may speculate that neutrinos are clustered 
around the stars and galaxies, or dark-matter clusters, in a pattern between photons 
and baryonic matter. In analogy to the air atmosphere surrounding the earth, there 
could be also neutrino atmosphere or clouds around stars and galaxies, or dark-matter 
clusters. As neutrinos interact with baryonic matter very weakly, the neutrino clouds 
can penetrate through stars and galaxies. As far as we know, there is no convincing 
measurement of neutrino clouds yet. It is thus necessary to design practical methods 
for the detection or constraint of the relic cosmic neutrino background (CVB), perhaps 
through detection of neutrino clouds of certain form. We will indicate in this work that 
systematic studies of high energy cosmic ray energy spectra of protons and nuclei may 
serve for such purpose. 

We first look at the energy spectrum of high energy cosmic ray protons, and 
check the influences due to interact with CMB photons and CVB neutrinos in their 
propagation. We know that there should be a big suppression in the spectrum, known 
as Greizen-Zatsepin-Kuzmin (GZK) cutoff jS], at threshold WL ~ 10 20 eV due to pion 
production by incident protons on CMB photons, p+7 — > n+N . Below that, there is also 
a suppression caused by electron-position pair production through p + 7 — > e + + e~ + p 
at threshold E^ e = 7.5 x 10 17 eV. These two processes will produce two "knees" due 
to suppressions in the proton energy spectrum above the corresponding thresholds, and 
one would also expect to see two "bumps" in the spectrum just below the thresholds 
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owing to the pileup of protons from higher energies. Compton scattering p + 7 — > 7 + p 
does not require an energy threshold, and the energy loss is also small for each collision, 
therefore it does not cause any non-smooth feature in the spectrum. We proceed to 
consider the interaction of the cosmic ray protons on the CuB neutrinos, 

p + v e — > e + + n. (1) 

When the proton with 4-momentum p = (E, p) interacts with the neutrino with 4- 
momentum k = (e u , k), and composes into a system with center of mass energy squared 
S, we have 

S — ml — ml , , 

E = 7 , P r, 2 

2e„ (l - ^l-ml/E 2 (\k\ /e„) costf) 

where 9 is the angle between p and k, m p is the proton mass, and m v is the neutrino 
mass. The energy of the proton E must be very large, so that we have 

S — m^ 

E ~ 2e„(l-(|k| /e^cosfl)' (3) 



The threshold energy is 

m 

Te„(l + |k| /e v ) 



m n + m e ) 2 - ml 



K e « vv, . ,^ p - ( 4 ) 



For zero mass neutrinos we have e v = (4/ll) 1 ^ 3 e 7 = 4.54 x 10~ 4 eV and |k| je v = 1, 
thus we get E^ e = 1.87 x 10 18 eV, which is between the pion and e + e~ production 
thresholds. Assuming the neutrino energy e u w m Ve ~ 1 eV, then |k|/e„ ~ 0, thus 
we get Ep e = 1.69 x 10 15 eV, which is well below the pion and e + e~ productions on 
CMB photons. Adopting the upper limit of neutrino mass m Ue < 2.2 eV from tritium 
experiments we get E^ e > 7.70 x 10 14 eV. Thus the threshold position E^ e ranges 
from 7.70 x 10 14 to 1.87 x 10 18 eV. In principle, the process ((H) should also produce a 
"knee" in the proton cosmic ray spectrum above the threshold and a "bump" below, as 
suggested by Wigmans [TH] . But the magnitudes might be very small, we need to check 
the constraints on CuB for any observable consequences. 

The mean energy of relic neutrinos is too small for interaction with the ordinary 
baryonic matter on earth, such that a high condensed neutrino clustering in cosmic 
background is not in conflict with our conventional knowledge. If relic neutrinos are 
condensed with a clustering factor c v compared to the mean number density, the number 
density should be n v = 56c u cm~ 3 . The cross section of reaction (JTJ) is of the order 

"~ 10 " 3 cm2 (^Hiw) 2 ~ ir42cm2 ' («) 

for m u ~ 0.5 eV and E v v = 4 x 10 15 eV calculated later on. Therefore the mean free 
path for a proton to propagate in Cz/B is 

1 4 x 10 15 

Ap E = — — ~ Mpc. (6) 

n,,a c,. 
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A clustering factor of c Ve ~ 10 13 would therefore correspond to a mean free path 
Ap e ~ 400 Mpc. As the protons with energies below the GZK cutoff and e + e~ production 
threshold can come from any source as large as the Hubble length, we thus predict an 
observable effect in the energy spectrum of proton cosmic rays, provided that relic 
neutrinos are clustered with a factor of the order 10 13 or larger. The cosmic neutrino 
clustering may occur around our galaxy, or around dark-matter clusters or supermassive 
black holes which are possible origin for the ultrahigh energy cosmic ray acceleration 
[Zj. The sizes of these neutrino clusterings are finite, might be a few orders larger than 
the sizes of the galaxies and dark-matter clusters. If we take the neutrino clustering 
size as the same order of a typical galaxy cluster size, i.e., 20 Mpc, then we will need a 
larger clustering factor of c v = 10 14 for the neutrino clustering to be observable. In this 
case, 30% of the primary cosmic ray protons and nuclei from outside the galaxy cluster 
would experience at least once interaction with cosmic neutrinos. 

We briefly review the available studies on the relic neutrino clustering. The 
clustering of relic neutrinos is suggested by Weiler in his proposal of using v + v — > Z 
to detect the relic neutrinos jH]. In order to explain the ultrahigh energy cosmic ray 
events above the GZK cutoff, the Z-bursts are suggested as the possible source, and the 
clustering factor for the cosmic neutrino background is adopted from a range of 10 2 ^ 7 
[HI El UD] to 10 8 ^ 14 [TT] to generate the required flux. A large neutrino clustering with 
order of 10 13 is also considered ^21 EI] to explain an anomaly in the tritium beta decay 
spectrum [T3*l IT4*] as from a background contribution of v e + 3 H — > e~ + 3 He. It is 
suggested by Wigmans ^H] that the knee of all-particle cosmic ray spectrum at energy 
4 x 10 15 eV is due to the process (JTJ) with a neutrino mass m Ve ~ 0.4 eV, provided 
with a dense neutrino clustering. There is also a theoretical calculation of the neutrino 
clouds as a result of very weakly interaction between neutrinos by the exchange of a 
very light scalar boson [T^j. Therefore a neutrino clustering factor in a range 10 2 ^ 14 
is by no means nonsense; it is still 16 — > 28 orders smaller than the baryonic matter 
clustering factor 4 x 10 30 . It is important to sharpen the range, and to reveal or to rule 
out the neutrino clustering with convincing evidence. We will show that the separation 
of the high energy cosmic rays between protons and nuclei can do this. 

Besides the pion and e + e~ productions as in the proton case, CMB photons can 
also disintegrate cosmic ray nuclei with energy of the order 10 19 eV and above through 
photonuclear processes such asA + 7^A^+(A — 1) |3|IT7j. There are similar processes 
for the interaction between Cz^B neutrinos and cosmic ray nuclei. Taking helium nuclei 
as an example, there are processes 



4 He + u e -> e+ + n + 3 H, 
4 He + v e -> e~ + p + 3 He, 



(7) 



which can be expressed in a general form as 
A + u e (u e ) ^ e ± + N + (A - 1). 



(8) 
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Thus we get the threshold energy for the cosmic ray nuclei 

E u e = S -m\ = (m A _i + m N + m e ) 2 - m\ 
A 2e„(l + |k| /e„) ~ 26,(1 + |k| /e v ) ' 1 ' 

If Cz/B neutrinos can cause an observable effect to the cosmic ray protons by p + u e — > 
e + + n at threshold energy Ep e , i.e., producing a "knee" above the threshold and a 
"bump" below in the cosmic ray proton spectrum, then these neutrinos should also 
produce a "knee" in the spectrum of cosmic ray A nuclei above E v £. There is no 
corresponding "bump" below the threshold, as the cosmic ray A nuclei with higher 
energies will be disintegrated into cosmic ray A-l nuclei or others. Comparing eqs. (@J) 
and (JIU), we get the relation between the threshold energies 

E v / = (m TO + m e f - m 2 p 

(m A _i + m N + m e ) 2 — m\ 

Thus systematic measurements of the spectra of cosmic rays for protons and nuclei can 
test whether process (0) is an observable effect in the cosmic rays or not. As there could 
be various kinds of nuclei in the cosmic rays, such as 12 C and 56 Fe, the corresponding 
thresholds can be exactly predicted with well established knowledge of nuclear physics. 
It should be able to measure or constrain the neutrino clustering factor c u , provided 
that the cosmic rays can be separated in composition between protons and nuclei. 

Once the threshold positions can be optimistically identified, we can also use these 
positions for a direct measurement of the neutrino mass, e.g., from eq. (@J) we get 

_ K + m e f - m 2 p 

Taking the starting point of the "knee" in all particle cosmic ray spectrum at energy 
4 x 10 15 eV as E p e , we get m Ve = 0.42 eV, which is within available experimental 
constraints jB] and is also consistent with the estimate by Wigmans |15j . One may 
naturally expect that protons are dominated in the measured all-particle cosmic rays 
around the knee, however, a claim of a so large effect of the knee as arising from Cz/B 
neutrinos ^H] seems extraordinary. The origin of the knee in the all-particle spectrum 
remains an open problem and there are also other novel explanations [TH] - It would be 
safe to check the corresponding knee in the energy spectrum for each species of cosmic 
ray nuclei. Assuming m Ue = 0.42 eV, we get the ratios of the threshold energies for p, 
4 He, 12 C, and 56 Fe cosmic rays, 

Up . £/4 He • £/12 C • -£/56 Fe / 12 x 

k 4x 10 15 : 1.8 x 10 17 : 4.4 x 10 17 : 1.3 x 10 18 . 1 ' 

Adopting a different value of m Ve will change these threshold positions, but do not 
change the ratios. Also the above positions should be considered as a rough estimate, 
as detailed consideration may change the position somewhat, for example, contribution 
from processes such as 

56 Fe + v e -»• e+ + 56 Mn, 

56 Fe + J/e ^ e - + 56 CO) ( 13 ) 
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where 56 Mn and 56 Co are unstable nuclei, may produce lower threshold positions for 
56 Fe at around 5 x 10 17 eV. 

Balloon and satellite investigations can provide direct measurements of the cosmic 
ray spectrum for each species of nuclei respectively. Unfortunately, the available 
measurements can only reach to the energy scale of 10 15 eV still 1 — >• 3 orders 
below the energy scale for our purpose. The available measurements of cosmic rays with 
energies ranging from 10 14 to a few 10 20 eV are constructed from extensive air showers 
[20] • The energy and composition of the primary cosmic rays are not directly measured, 
but can be reconstructed from the air shower profiles such as electron and muon size 
distributions [TBI EH - 

There have been experimental attempts to de-convolute the all-particle energy 
spectrum around the knee into four individual primary mass groups of p, He, C, and 
Fe by the KASCADE collaboration |18j . The mass information of primary extensive air 
showers is deduced by the secondary electron and muon shower size distributions. Each 
of these preliminary energy distributions of the four mass groups exhibits a knee like 
structure with an increase of their knee positions for heavier A nuclei. It is interesting to 
note that the observed positions are in agreement with the above predictions, as can be 
checked from the KASCADE results ^H]- It is also possible that the proton knee starts 
at (1 — > 2) x 10 15 eV [221, which corresponds to the neutrino mass m Ue = 0.84 — > 1.68 eV, 
so that all corresponding knee positions from Eq. ()12j) are shifted by a factor 2 — > 4 to 
the lower energy side. Considering that there are large uncertainties in both theoretical 
predictions and experimental reconstructions, we may take the KASCADE observation 
as possible support for a large neutrino clustering factor to be measurable. Of course, 
more precision and reliable measurements of detailed energy spectra of cosmic ray 
protons and nuclei are still necessary, before one can draw a conclusion to confirm 
or rule out such a situation. 

The relative abundances of cosmic ray nuclei may show more significant signature 
than the isolated spectrum for each species of nuclei. Thus we may have chance to 
detect and constrain the CuB with a smaller clustering factor. For example, the ratio 
of 3 He/ 4 He is of the order 10~ 4 in the universe from primordial nucleosynthesis j2H|. If 
this feature is also true at the source for very high energy cosmic rays, then interaction 
of abundant cosmic ray 4 He nuclei with cosmic neutrinos, i.e., 4 He + v e — *> e~ + p + 3 He, 
will produce secondary cosmic ray 3 He nuclei. The energy of the produced 3 He can 
be exactly calculated and it should be below and around the 4 He threshold energy. 
This will cause a significant "bump" in the ratio of 3 He/ 4 He around the corresponding 
energy position and a "tail" below, in case the neutrino clustering factor is of the order 
10 9 or larger. Contributions from CMB photon disintegration of abundant cosmic ray 
nuclei should happen at higher energy scale, thus do not diminish the signature due to 
interaction of the cosmic ray nuclei with cosmic neutrinos. One may also search for a 
pair of nuclei with more significant difference in their relative abundances. This can 
reduce the observable neutrino clustering factor to even much smaller cases. 

The abundance of cosmic ray nuclei can show more significant effect in case with a 
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higher neutrino clustering factor. For example, the number of less abundant nuclei will 
be enhanced below and around the knee as secondaries resulting from collisions of the 
primary abundant cosmic ray nuclei with CzaB neutrinos. We can use this aspect as a 
further signature for neutrino clustering, in addition to the knee onset position for each 
species of cosmic ray nuclei. Other models for the knee may also produce an increase of 
knee position for cosmic ray nuclei with larger A ^H], and it is rather difficult to make a 
distinction for the knee positions with charge number Z dependence or nuclear number 
A dependence [21]. But the energy dependence in the relative abundances of cosmic ray 
nuclei should be different in different models. Thus the cosmic ray spectra of protons and 
nuclei with detailed energy dependence and abundance information can provide decisive 
detection or constraint on the cosmic neutrinos. Indeed, the isotopic ratio of 3 He/ 4 He has 
been measured in balloon experiments (2E], and this ratio is found significantly enhanced 
compared to that from nucleosynthesis, with also a trend to increase as energy increases. 
However, the energies of those cosmic rays only range around 10 8 to 10 10 eV, rather low 
to conclude as the "tail" resulting from the mechanism we suggested. 

In comparison with Wigmans' proposal to attribute the knee in high energy cosmic 
ray spectrum as arising from CzaB neutrinos, we provide a practical relation, i.e., 
Eq. (|12jh between the "knees" caused by interactions of Cz/B neutrinos with the protons 
and specific species of nuclei ||. Such relations can provide a crucial test of the idea to 
attribute the existing knee as a signal for the interaction between the protons with Cz/B 
neutrinos. We also propose a practical method to measure the isotopic ratio such as 
3 He/ 4 He, so that we may have chance to detect and constrain the Cz/B with a smaller 
clustering factor. 

In summary, we proposed a method to investigate the scenario that cosmic relic 
neutrinos are highly clustered around stars and galaxies, or dark-matter clusters, rather 
than uniformly distributed in the universe. High energy cosmic ray protons and nuclei 
should interact with these cosmic neutrinos through their propagation in the universe. 
There should be observable signature in cosmic ray energy spectra of protons and nuclei, 
for a neutrino clustering factor beyond 10 13 . The onset positions of the knee signature 
should be closely related between protons and nuclei, and there seems some support 
from existing experiments. We have also indicated that a smaller neutrino clustering 
factor can be detected or constrained by an increase of secondary less abundant cosmic 
ray nuclei produced from interaction of the corresponding primary abundant cosmic ray 
nuclei with the cosmic neutrino background. The neutrino mass can also be directly 
measured from the knee onset positions reflecting interactions of the high energy cosmic 
ray protons and nuclei with cosmic neutrinos. Thus systematic investigations of high 
energy cosmic ray protons and nuclei are practical to reveal or constrain the cosmic 
neutrino clustering and neutrino properties. Precision measurements are required for 
this purpose in balloon, satellite, and extensive air shower experiments of high energy 
cosmic rays with explicit nuclei composition information included. 

|| The discussion on the reaction of CVB neutrinos with a particles by Wigmans is not appropriate 
from nuclear physics, as the corresponding reaction should be Eq.(7) in this paper 
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